STRUCTURED THALAMOCORTICAL CONNECTIVITY REVEALED BY 
RANDOM WALKS ON COMPLEX NETWORKS 
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The segregated regions of the mammalian cerebral cortex and thalamus form an extensive and 
complex network, whose structure and function are still only incompletely understood. The present 
article describes an application of the concepts of complex networks and random walks that allows 
the identification of non-random, highly structured features of thalamocortical connections, and their 
potential effects on dynamic interactions between cortical areas in the cat brain. Utilizing large-scale 
anatomical data sets of this thalamocortical system, we investigate uniform random walks in such a 
network by considering the steady state eigenvector of the respective stochastic matrix. It is shown 
that thalamocortical connections are organized in such a way as to guarantee strong correlation 
between the outdegree and occupancy rate (a stochastic measure potentially related to activation) 
of each cortical area. Possible organizational principles underlying this effect are identified and 
discussed. 



The relationship between the topology of structural 
connection patterns and cortical dynamics currently rep- 
resents a significant challenge to brain theory. Numerous 
neuroanatomical studies have revealed that the pathways 
of the mammalian thalamocortical system exhibit spe- 
cific patterns ranging in scale from interconnections link- 
ing whole brain regions to intra-areal patterns of connec- 
tions between cell populations or individual cortical neu- 
rons. Comprehensive descriptions of large-scale anatom- 
ical patterns of cortical connectivity have been collated 
for several mammalian species (e.g. 0,0)- Quantitative 
analysis has revealed that these patterns are neither com- 
pletely regular nor completely random but exhibit 
specific attributes, such as short path lengths combined 
with high clustering 13, 0] , short wiring and distinctive 
hierarchical features [(j. 

The current article proposes a novel macroscopic ap- 
proach to studying neuronal dynamics in cortical maps 
based on concepts from random walks (e.g. 0, 111) an( i 
complex networks (e.g. [UE1E1])- Specifically, we ex- 
amine the interactions between cortex and thalamus in 
the mammalian central nervous system. The cerebral 
cortex consists of a network of interconnected function- 
ally specialized regions, each of which maintains recip- 
rocal connections with a specific set of thalamic nuclei. 
These nuclei relay cortical inputs and outputs, and are 
thought to play a role in coordinating cortical interac- 
tions. The thalamocortical architecture is represented in 
terms of a complex network, obtained by assigning a node 
to each cortical region and thalamic nucleus (nodes) and 
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distributing anatomical connections as links (or edges) 
between such nodes. The flow of activation of each node 
is then modeled by an agent engaged in a random walk 
on the edges of the network. The movement of the agent 
expresses the exchange or flow of information between 
adjacent nodes, analogous to a diffusive process. After a 
long number of time steps T, the occupancy rate, i.e. the 
number of times each node was visited divided by T, cap- 
tures the rate of participation of each node in network- 
wide information flow, a dynamic measure analogous to a 
cortical activation level. In order to relate this dynamical 
measure to the underlying network topology, we examine 
the outdegree of each network node. 

The anatomical connection matrix of cat cortex and 
thalamus was obtained from the study by Scannell et 
al. |2| and contained 53 cortical regions and 42 thalamic 
regions (Fig. la). Most corticocortical connections and 
all thalamocortical connections are reciprocal (i.e. the 
matrix is nearly symmetric), and no connections exist 
between any of the thalamic nodes. Following the anal- 
ysis in Q the 53 cortical areas can be divided into four 
functionally distinct clusters: visual (17, 18, 19, PLLS, 
PMLS, ALLS, AMLS, VLS, DLS, 21a, 21b, 20a, 20b, 7, 
AES, PS), auditory (AI, All, AAF, P, VPc, EPp, Tern), 
somatomotor (31, 3b, 1, 2, SII, SIV, 4g, 4, 61, 6m, 5Am, 
5A1, 5Bm, 5B1, SSAi, SSAo), and frontolimbic (PFCMil, 
PFCMd, PFCL, la, Ig, CGa, CGp, RS, 35, 36, pSb, Sb, 
Enr, Hipp). Visual and auditory clusters are combined 
into a 'posterior' cluster, while somatomotor and fron- 
tolimbic clusters are combined into an 'anterior' cluster. 
The adjacency matrix is henceforth represented as K, 
with the presence of a directed connection from region j 
to i indicated as K(i,j) = 1. 

Because of the relatively small number of involved re- 
gions and the fact that any node can be reached from 
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any other node in the network (i.e. it forms a single 
connected component), it is possible to characterize the 
steady state node occupation effectively in terms of the 
dominant eigenvector of the respective stochastic matrix 
(the random walk can be understood as a random driven 
Markov chain) I n order to do so, a stochastic version 
S of the adjacency matrix K is obtained by making 

0{j)=Y J K(i,j) (1) 

S(i,j)=K(i,j)/0(j) (2) 

where O(j) is the outdegree of node j, i.e. the number 
of edges emanating from that node. Note that 5 is a 
stochastic matrix because S(i,j) > for any i and j 
and J2iLiS(i,j) = 1 for any j. At steady state, the 
occupation rate of the random walk movements can be 
obtained from the eigenvector equation Sv — v, where v 
is the eigenvector of the stochastic matrix S associated to 
the unit eigenvalue . The eigenvector v is normalized 
so that 52»=i Vi = 1. The average occupancy rate for 
each cortical region i is therefore given by the eigenvector 
element v(i). 

Figure lb shows the outdegree (O) versus the occu- 
pancy rate (R) for the cortical regions obtained by the 
eigenvalue analysis considering all thalamocortical inter- 
connections. Posterior and anterior regions are identi- 
fied by 'x' and '*' respectively. This network exhibits a 
pronounced positive correlation (Pearson correlation r = 
0.84) between the topological and dynamical measures. 
We examined the hypothesis that this strong correlation 
is due to the specific organizational pattern of thalamo- 
cortical connections. First, we repeat the analysis con- 
sidering only the cortico-cortical connections. Now, a 
markedly weaker correlation (r = 0.54) was observed be- 
tween the topological (i.e. outdegree) and dynamical (i.e. 
occupancy rate) measures, which are illustrated in Figure 
lc. When replacing the thalamocortical connections by 
the same number of randomly distributed edges (main- 
taining reciprocity), we again obtained a weaker correla- 
tion (r = 0.65±0.044 considering 20 random realizations) 
between outdegree and occupancy ratio. Randomization 
of the corticocortical connections while keeping the orig- 
inal thalamocortical edges results in a decrease, though 
slightly less marked, of the correlation between outdegree 
and occupancy rates (r — 0.67 ± 0.041 for 20 random re- 
alizations). Similar figures were obtained regarding the 
correlation between indegree and activation rate. Such 
results seem to indicate that the specific pattern of tha- 
lamocortical connections promotes a strong correlation 
between topological and dynamical features of cortical 
interactions. 

In order to try to infer which structural features of the 
thalamocortical interconnections are principally respon- 
sible for generating the correlation between topology and 
dynamics, each thalamic node was isolated together with 
the respective in and outbound connections - defining 



a V-structure, as allowed by the complete lack of inter- 
connectivity between these nodes. The shortest paths 
between each of the cortical regions at the destination 
of the outbound edges of thalamic node i and all the 
cortical regions at the origin of the inbound edges are 
computed. The respective average (Av) and standard 
deviation values (St) are calculated and shown in Figure 
Id in decreasing order. The plot shows that the addi- 
tional edges (and therefore outdegree) imposed by the 
thalamic V-structures tend to implement cycles of just a 
few edges among the group of cortical regions to which it 
is attached. In other words, if the shortest path between 
two cortical areas attached to a V-structure involves k 
edges, the connection to the thalamic structure will de- 
fine a cycle with k + 2 edges. This effect contributes to 
a directly proportional increase of the occupancy rate, 
being potentially related to the observed correlation be- 
tween the considered topological and dynamical features. 
Other elements of the thalamocortical structure also con- 
tributing to the strong correlation include the marked 
similarity between the in and outdegree at each node as 
well as the fact that the whole network defines a con- 
nected component. 



In conclusion, we have shown in this article that a novel 
combination of concepts from complex networks and ran- 
dom walks/Markov chains can reveal important struc- 
tural properties underlying thalamocortical architecture 
in the cat brain. Our analysis demonstrates the presence 
of a strong correlation between network topology (in and 
outdegree) and dynamic features of cortical activation 
(occupancy rate) as modeled by random walks. The ba- 
sic feature of the thalamic connections that is responsible 
for establishing such a correlation has also been investi- 
gated in terms of shortest path analysis between cortical 
nodes participating in each of the thalamic V-structures. 
A possible critical feature of such circuits was identified 
as the reliatively short paths between subgroups of corti- 
cal regions. In a naive analogy, it is as if the thalamic V- 
structures provided a feedback mirror to groups of closely 
connected cortical regions. As for the possible implica- 
tions of the correlation between outdegree and occupancy 
rate, it has the two following interesting effects: (1) hubs 
of connectivity are likely to be hubs of neuronal informa- 
tion flow and activity; and (2) in the case of scale free 
topology (degrees), the dynamic participation or cortical 
activation at each node will also become scale free. 



The reported findings have potentially important im- 
plications for the understanding of cortical functional ar- 
chitecture and pave the way to a number of subsequent 
studies, including the comparison of cortical connectiv- 
ity of other species and the use of other kinds of random 
walks and dynamics. At the same time, several other 
problems in complex networks can be addressed in terms 
of the correlation measurement proposed in this article. 
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Adjacency matrix of the directed complex network de- 
fined by the cat thalamocortical connections (a). The 
scatterplot of outdegree (0) and occupancy rate (i?) con- 
sidering all thalamocortical connections (b) and only cor- 
tical connections (c). The average ± standard deviation 
of the shortest paths between the cortical nodes connect- 
ing to thalamic nodes, considering only cortical intercon- 
nections (d). 
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FIG. 1: Figure 1. 



